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1.  Introduction 

While  both  direct  numerical  simulation  (DNS)  and  large  eddy  simulation 
(LES)  are  at  or  near  the  top  of  a hierarchy  of  solution  methods  for  tur- 
bulent flow  fields,  the  Reynolds  number  constraints  and  extensive  compu- 
tational requirements  preclude  their  widespread  use  as  a practical  tool  in 
aerodynamic  applications.  Nevertheless,  these  methodologies  can  be  used 
in  support  of  more  practical  engineering  tools  such  as  Reynolds-averaged 
Navicr-Stokes  (RANS)  formulations. 

Aerodynamic  flow  fields  present  several  challenges  even  for  the  most  ro- 
bust RANS  approaches.  Geometric  complexity,  Mach  number  effects,  and 
Reynolds  numbers  of  0(1O7)  contribute  significantly  to  the  stringent  re- 
quirements needed  for  successful  flow  field  computations.  Nevertheless,  even 
with  this  apparent  disparity  between  the  capabilities  of  DNS  and  LES  and 
the  requirements  for  accurate  prediction  of  relevant  aerodynamic  flow  fields, 
both  types  of  simulations  can  provide  useful  information  if  properly  chosen 
“unit  problems”  are  studied.  Two  examples  of  practical  aerodynamic  flow 
fields  will  serve  to  highlight  some  of  the  critical  dynamical  problems  as  well 
as  the  need  for  well  chosen  simulations  that  can  help  improve  the  predictive 
capability  of  RANS  calculations  or  composite  solution  approaches. 

The  inability  of  RANS  formulations  to  correctly  predict  critical  dynamic 
features  of  complex  flow  fields  lies  with  inadequacy  of  the  closure  models 
for  the  higher-order  correlations  that  appear  in  the  RANS  formulation.  A 
detailed  discussion  of  such  closures  is  outside  the  scope  of  the  current  topic; 
however,  correlations  involving  the  pressure-strain  rate  or  pressure  velocity, 
triple-velocity  and  (tensor)  dissipation  rate  all  can  contribute  significantly 
to  the  predictive  accuracy  of  a RANS  formulation.  Since  experimental  stud- 
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ies  are  unable,  in  general,  to  accurately  measure  these  higher-order  correla- 
tions, numerical  simulations  that  can  help  delineate  the  role  of  such  terms 
in  complex  flows  is  extremely  useful. 

Ever  increasing  computational  power  as  well  as  improved  understand- 
ing of  the  effect  of  filter  cut-off'  on  subgrid  scale  models  will  continue  to 
move  LES  (and  DNS)  towards  the  realm  of  an  engineering  tool.  In  ad- 
dition, extensions  of  RANS  formulations  to  unsteady  RANS  and  related 
methodologies  will  continue  to  move  these  type  of  formulations  closer  to- 
ward the  realm  of  a large-eddy  solution  technique.  The  future  of  obtaining 
such  composite  formulations  will  depend  on  the  inherent  consistency  be- 
tween the  (velocity)  fields  computed.  An  example  of  such  a consistency 
requirement  will  be  discussed  here  and  its  effect  on  model  development. 

2.  Airframe  Dynamics 

High- lift  devices  composed  of  multi-element  airfoils  (Fig.  1)  generate  com- 
plex aerodynamic  flow  fields.  At  cruise  conditions,  these  flows  are  generally 
at  the  high  subsonic  speed  range  although  conditions  can  exist  in  some  re- 
gions for  shocks  to  develop.  Each  element  of  the  multi-element  system  has 
characteristic  features  which  are  of  interest  in  their  own  right  (Rumsey  et 
al. , 1998). 
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Figure  1.  Sketch  of  multi-element  airfoil  configuration  and  important  dynamic  charac- 
teristics (C.  L.  Rumsey,  private  communication) 


Common  to  all  elements  is  the  need  for  accurate  (natural)  transition  pre- 
diction. Within  a RANS  formulation,  the  turbulence  closure  models  need 
to  be  properly  sensitized  to  the  transition  process.  Corresponding  to  the 
usual  practice  of  turbulence  model  calibration  using  both  DNS  and  LES 
results,  can  such  simulations  be  performed  on  (natural)  transitioning  flows 
with  the  goal  of  providing  the  information  needed  for  calibration  of  RANS 
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closure  models?  Simulations  of  transitional  flows  have  been  undertaken  pre- 
viously; however,  the  motivation  was  not  to  develop  a data  base  for  model 
development  but  was  to  obtain  a detailed  description  of  the  (nonlinear)  sta- 
bility process  leading  to  a turbulent  flow.  As  such,  parameter  ranges  and 
dynamics  studied  did  not  focus  on  the  data  needed  to  help  refine  models  for 
disturbance  stress  fields  which  would  be  suitable  for  RANS-type  closures. 

Another  feature  of  this  multi-element  system  is  that  all  the  elements  are 
dynamically  coupled  through  the  interaction  of  a wake(s)  generated  by  an 
upstream  element  and  the  boundary  layer  of  a downstream  element.  It  has 
not  been  possible  to  accurately  predict  the  downstream  evolution  of  these 
generated  wakes  and  as  such,  the  downstream  predictions  have  not  been 
accurate  as  well  (Rumsey  and  Gatski,  2002).  It  is  not  clear  whether  the  clo- 
sure model  deficiencies  are  due  to  the  strong  pressure  gradients  imposed  on 
these  wakes  (such  as  the  slot  generated  wake),  or  the  sequence  of  pressure- 
gradients  the  developing  wakes  experience  as  they  evolve  downstream.  In 
isolation,  prediction  of  the  near-wake  flow  field  has  been  successful,  at  least 
over  a limited  range  of  pressure  gradients  (Carlson  et  al.,  2000).  While 
DNS  of  a wake  flow  has  been  performed  previously  (Moser  et  al.,  1998), 
it  did  not  focus  on  effects  of  pressure  gradient  or  other  issues  related  to 
the  multi-element  dynamics  of  interest  here.  Nevertheless,  a well-posed  and 
accurate  DNS  (or  LES)  could  be  used  to  validate  the  predictive  capability 
of  a RANS  model  for  such  wakes  in  relative  isolation  to  other  complicating 
dynamic  features  associated  with  the  full  configuration. 

In  addition  to  transition  location  and  wake  evolution  which  affect  the 
global  behavior  of  the  multi-element  flow  field,  each  element  of  the  system 
also  has  some  characterizing  dynamic  features.  The  slat,  for  example,  is 
composed  of  a curved,  transitional  flow  which  can  sustain  shocks  depend- 
ing on  configuration,  and  a possible  unsteady,  transitioning  shear  flow  in  the 
lower  surface  cove  region.  Thus,  shock  boundary-layer  interaction  predic- 
tion on  the  curved  upper  surface  can  be  an  important  feature  of  any  overall 
prediction  scheme  of  high-lift  dynamics.  Unfortunately,  the  ability  to  accu- 
rately predict  such  shock-boundary  layer  interactions  are  strongly  problem 
dependent.  On  the  lower  surface  of  the  slat  experiments  have  shown  an  in- 
herent unsteadiness  in  the  flow.  Unfortunately,  such  measurements  are  dif- 
ficult and  detailed  mappings  are  unavailable.  In  any  case,  such  unsteady  ef- 
fects produce  a significant  challenge  to  any  prediction  scheme.  Well-focused 
simulations  of  model  flows  which  can  partially  capture  some  of  these  key 
features  of  the  practical  configuration  can  be  critical  to  identifying  specific 
deficiencies  of  a particular  RANS  closure  scheme. 

The  main  element  has  been  less  problematic  once  the  transition  location 
is  properly  predicted  or  fixed.  Errors  in  the  slat  wake  predictions  inherently 
persist  downstream;  however,  the  slat  wake  velocity  deficit  only  begins  to 
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interact  with  the  main  element  boundary  layer  near  the  aft  portion  of  the 
element.  In  addition,  the  flow  over  the  flap  is  moderately  curved  which 
further  complicates  the  dynamics.  Turbulence  models  currently  available 
are  only  moderately  successful  in  capturing  all  these  dynamics,  even  in  well- 
posed  comparative  studies  with  experiments.  Clearly,  a significant  challenge 
is  to  identify  suitable  test  flow  fields  that  can  be  accurately  computed  with 
DNS  or  LES  approaches,  and  still  provide  a database  that  can  be  used  to 
validate  turbulent  closure  models  that  could  lead  to  improved  closures. 

The  mean  velocity  field  over  the  flap  is  characterized  by  a velocity  deficit 
produced  by  the  merging  of  the  slat  wake  with  the  main  element  boundary- 
layer  flow  and  an  altered  inner  layer  velocity  field  due  to  the  gap  between 
the  main  element  and  flap.  As  the  angle-of-attack  increases,  flow  separation 
over  the  high-lift  device  may  be  initiated  over  the  flap  adding  to  the  com- 
plexity of  the  flow.  Obviously,  prediction  schemes  are  not  successful  over  a 
wide  parameter  range  and  could  benefit  from  simulations  which  isolate  key 
features.  One  area  of  particular  interest  would  be  the  flow  field  dynamics 
in  the  region  downstream  of  the  separation  point.  A detailed  mapping  of 
this  flow  would  hopefully  allow  for  improved  model  development. 

3.  Engine  Flow  Fields 

A second  example  of  a complex  aerodynamic  flow  field  is  that  associated 
with  scramjet /ramjet  engines  (Fig.  2).  The  inlet  flow  field  is  composed 
of  complex  shock  boundary-layer  interactions  and  the  detailed  dynamics  of 
such  flows  is  not  well  understood.  Downstream  in  the  combustor,  high  speed 
fuel  injection  occurs  resulting  in  a complex  mixing  process  with  shocks 
present. 
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Figure  2.  Sketch  of  inlet  flow  field  of  scramjet/ramjet  engine  and  important  dynamic 
characteristics  (J.  A.  White,  private  communication) 

Once  again,  such  flow  field  predictions  are  beyond  the  capabilities  of  ei- 
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ther  direct  or  large-eddy  simulations.  Thus,  at  this  time  a RANS  approach 
is  the  only  practical  means  of  obtaining  predictions  of  the  flow  field.  The 
challenge,  then,  is  to  identify  important  dynamical  features  which  charac- 
terize different  regions  of  the  flow  field  and  formulate  suitable  DNS  and 
LES  unit  problems  that  can  be  studied  to  improve  the  RANS  closure  mod- 
els. Compressible  RANS  approaches  using  Favre-averaged  variables  have 
been  the  methodology  most  commonly  used  in  such  flows.  Such  approaches 
yield  additional  correlations  in  the  governing  equations  which  are  due  to 
compressibility  and  for  which  closure  models  need  to  be  developed.  These 
may  include  heat  and  mass  flux  models  as  well  as  dilatation  models.  De- 
pending on  the  degree  of  compressibility,  the  turbulence  can  be  significantly 
altered  from  its  incompressible  behavior.  If  so,  several  modeling  issues  arise 
which,  unlike  the  incompressible  case,  need  to  be  resolved.  While  numerical 
simulations  have  been  performed  in  the  last  decade  on  compressible  flow 
fields  (e.g.  Freund  et  al.  2000,  Sarkar  and  Pantano  1999),  the  results  have 
not  been  extensively  utilized  in  model  development.  Thus,  at  the  outset  it 
would  be  beneficial  to  examine  simulations  already  performed  and  evaluate 
their  usefulness  on  improving  existing  currently  available  models. 

As  Fig.  2 shows,  the  inlet  region  is  dominated  by  a complex  shock  pat- 
tern which  has  a significant  impact  on  the  flow  field.  Useful  information 
could  be  obtained  from  a simple  boundary-layer  computation  with  an  im- 
pinging shock  (e.g.  Adams  2000).  Such  a calculation  would  give  some  insight 
into  the  dominant  dynamic  features  that  would  need  to  be  considered  in  the 
full  problem.  Another  important  feature  of  such  flows  is  the  high-speed  fuel 
injection  process.  The  injection  produces  a very  complex  flow  field  which  is 
very  challenging  to  predict.  Can  a representative  DNS  or  LES  unit  problem 
be  devised  that  would  successfully  replicate  some  of  the  important  physics? 
Such  a database  could  be  used  to  validate  existing  models  or  lead  to  the 
development  of  improved  closure  models. 

Certainly  complex  flow  field  dynamics  are  not  limited  to  the  scram- 
jet/ramjet  engines.  Unlike  the  scramjet/ramjet  engines,  turbofan  engines 
contain  a rotor/stator  assembly  which  introduces  extra  strain  effects  into 
the  turbulent  flow  field.  Flow  curvature  and  non-inertial  effects  adds  to  an 
already  complex  flow  that  includes  leading-edge  shocks  ahead  of  the  rotor 
blades,  shedding  vortices,  and  unsteady  wakes. 

4.  Future  Challenges 

Many  new  challenging  problems  continue  to  arise  predicated  on  the  need  to 
better  control  the  flow  field  dynamics  in  complex  configurations.  Adaptive 
flow  and  noise  control  are  but  two  examples  along  with  the  continuing  need 
to  improve  the  predictions  of  the  type  of  flow  fields  discussed  in  the  previous 
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sections. 

Adaptive  or  active  flow  control  comprise  a closed  system  where  a con- 
tinuous feedback  loop  exists  to  optimize  some  flow  field  characteristic.  Ac- 
tive cavity,  fluid/structural  shaping  and  separation  are  examples  of  control 
mechanisms  that  can  be  used.  Each  of  these  control  mechanisms  when 
embedded  in  complex  aerodynamic  configurations  would  require  highly  re- 
solved simulations  in  order  to  adequately  describe  the  complex  flow  field 
dynamics.  The  same  holds  true  for  improved  airframe,  fan  and  jet  noise 
control.  Detailed  simulations  of  full  configuration  slat  leading  edge,  flap 
trailing  and  side-edge,  and/or  landing  gear  geometries  are  prohibitive. 

As  in  the  previous  sections,  DNS  and  LES  methodologies  may  be  best 
suited  to  problems  that  describe  the  complex  flow  field  structure  in  simpli- 
fied geometries  but  which  still  capture  the  key  dynamic  features  of  the  full 
flow  field.  Another  approach  would  be  to  develop  composite  methodologies 
capable  of  achieving  the  accuracy  of  the  DNS  and  LES  methodologies,  but 
in  full  configurations.  Such  composite  or  hybrid  approaches  have  begun  to 
be  formulated  and  these  include  the  Detached  Eddy  Simulation  (DES)  (e.g. 
Spalart  1999)  and  the  Flow  Simulation  Methodology  (FSM)  (e.g.  Zhang  et 
al.  2000). 

Such  composite  approaches  can  be  constructed  from  modified  Reynolds 
averaged  Navier-Stokes  (RANS)  and  LES  formulations,  for  example.  The 
usual  RANS-type  formulations  are  probably  unacceptable  since  they  do 
not  in  general  handle  non-equilibrium  effects  properly  and  are  not  properly 
sensitized  to  the  broad  spectrum  of  scales  present.  Extensions  to  the  usual 
RANS-type  formulations  include,  for  example,  the  triple-decomposition  ap- 
proach originally  proposed  by  Reynolds  and  Hussain  (1972)  and  the  Semi- 
Deterministic  Method  (SDM)  (Ha  Minh  and  Kourta  1993).  Such  extensions 
inherently  solve  time-dependent  RANS-type  formulations  and  as  such  have 
been  referred  to  as  unsteady  RANS  (URANS),  time-dependent  RANS  (T- 
RANS)  and  VLES.  The  label  VLES  (Very  Large  Eddy  Simulation)  of  course 
being  put  forward  to  establish  a formal  link  with  the  LES  approach.  Some 
new  alternatives  to  LES  are  now  also  appearing  such  as  the  Navier-Stokes- 
o.  model  which  includes  nonlinear  dispersive  effects  (Chen  et  al.  2000),  and 
the  Coherent  Vortex  Simulation  (CVS)  approach  (Farge  et  al.  2000,  see  also 
Goldstein  et  al.  2000).  It  remains  to  be  seen  whether  such  approaches  pro- 
vide better  modeling  insight  for  RANS-type  closures  or  a more  conducive 
basis  for  composite  formulations. 

It  should  be  recognized  that  while  the  RANS  and  LES  equations  are 
formally  equivalent,  the  flow  field  is  being  described  differently  due  to  the 
disparity  of  scales  being  resolved.  Thus,  one  might  naturally  ask  whether 
there  is  some  way  in  which  both  methods  would  produce  the  same  flow  field 
as  described  by  the  velocity  field,  for  example. 
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It  is  worthwhile  to  look  at  a strategy  proposed  by  Germano  (1999). 
Consider  the  ensemble  mean  of  both  the  velocity  and  pressure  fields  and 
assume  that  for  /(x,  t)  (?q(x,  i)  or  p(x,i)) 


#{/(M)}  = £{/(M)}  = £(M) 

DNS  LES  ANS 


(1) 


The  ambiguity  associated  with  averaged  Navier-Stokes  variable  is  inten- 
tional. In  both  DNS  and  LES,  the  statistical  average  of  the  variable  (DNS) 
or  the  filtered  variable  (LES)  is  performed  simultaneously  with  (or  post- 
processed  after)  the  numerical  solution  of  the  Navier-Stokes  equations.  In 
RANS-type  formulations  (which  are  labeled  here  as  ANS  for  Averaged 
Navier-Stokes),  it  is  the  averaged  variable  which  is  solved  for  directly.  For 
this  reason,  the  ANS  mean  velocity  and  pressure  will  simply  be  represented 
by  Ui  and  P,  respectively. 

Now,  in  addition  to  Eq.  (1)  which  assumes  the  equivalence  of  the  en- 
semble mean  of  the  filtered  quantity  with  the  ensemble  mean,  it  is  also 
required  that  the  corresponding  mean  momentum  equations  be  equivalent, 
then 
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Note  that  here  the  second-moment  of  the  fluctuating  velocity  rt]  is  given 
by  some  statistical  average  of  (ul  — Ui ) (uj  — Uj). 

The  partitioning  suggested  by  <Sles  can  also  be  clearly  shown  in  spec- 
tral space.  For  a statistically  homogeneous  field  using  a sharp  cut-off  filter 
(G(k)  = H(kc-\k\)) 

E{(ui  -E{ui})  (uj  -£{«,})}  = f d3 k G2(k)$„(M) 

J V/c 

= f d3k$y(M)  (9) 
J\k\<kc 

and 

E{UiUj  - UiUj}  = E{rffs}  = f d3k  (1  - G(k))2$y  (k,  t) 

JVk 

= f d3k$ij{k,t)  (10) 

where  the  energy  spectrum  tensor  (k,  t)  is  related  to  the  Fourier  trans- 
formed instantaneous  velocity  U{  by 

E{(ui  - E{ui })  (uj  - E{u-j})}  = 8\ k + k')$jj(k,  t)  (11) 

Equations  (9)  and  (10)  show  that  as  the  cut-off  wavenumber  increases,  the 
contribution  from  r^GS  diminishes  and  the  LES  formulation  evolves  toward 
a full  DNS;  whereas,  as  the  cut-off  wavenumber  decreases,  the  contribution 
from  rffs  increases  and  more  of  spectrum  needs  to  be  modeled. 

In  order  that  the  velocity  fields  E{Tii}  and  Ut  computed  from  either  an 
LES  Eq.  (3)  or  a RANS-type  Eq.  (4)  formulation  be  the  same,  a formal 
requirement  for  consistency  between  the  two  methodologies  would  be 

nj  = E{ (ui  -£{”,})  (uj  -£{“j}) } + E{t f/s},  (12) 

or  using  Eqs.  (9)  and  (10), 

/ d3k$f,(k,i)  = [ d3k$i,(k,t)  + [ d3k$tj{k,t).  (13) 

J'ik  j|fc|<fcc  ^|fe|>fcc 

The  relation  in  (12)  can  also  be  shown  to  be  a direct  consequence  (Germano 
1999)  of  the  initial  assumption  about  the  mean  fields  in  Eq.  (1).  Germano 
(1999)  used  the  result  in  Eq.  (12)  to  derive  a simple  relation  for  the  subgrid 
scale  eddy  viscosity.  An  alternative  view  would  be  to  use  (12)  as  a guide  to 
improved  modeling  of  Tly 

In  a statistically  steady  flow,  where  the  ANS  is  simply  the  usual  RANS 
(e.g.  long  time  average)  formulation,  contains  the  entire  effect  of  the 


DNS/LES  FOR  NASA  AERODYNAMIC  NEEDS 


33 


turbulence  on  the  mean  flow,  and  in  a sense  corresponds  to  the  limiting 
value  for  a decreasing  cut-off  wavenumber  kc.  However,  in  statistically  un- 
steady flows  where  methods  such  as  unsteady  RANS  (URANS  or  T-RANS, 
VLES)  may  be  needed,  contributions  from  both  Eqs.  (9)  and  (10)  need  to 
be  properly  represented  by  in  order  for  consistency  with  an  LES  for- 
mulation. In  practice,  for  inhomogeneous  flows  the  mean  (velocity)  fields 
obtained  from  a RANS-type  formulation  are  associated  with  temporal  av- 
erages. This  further  complicates  any  consistency  arguments  due  to  the  fact 
that  any  coupling  between  such  a temporal  average  and  the  spatial  filtering 
associated  with  an  LES  would  necessarily  involve  a complicated  (and  prob- 
ably unknown)  dispersion  relation  (Pruett,  2000)  that  needs  to  be  taken 
into  account. 

5.  Summary 

While  direct  and/or  large  eddy  simulations  are  probably  not  going  to  be  ca- 
pable of  predicting  complex  aerodynamic  flow  field  themselves,  the  method- 
ologies can  and  should  be  utilized  to  solve  flow  problems  that  replicate  the 
essential  dynamic  features  of  the  full  problem.  In  the  near  term,  this  may  be 
the  most  useful  role  of  such  methodologies  rather  than  as  a substitute  to  the 
currently  available  RANS-type  models.  These  simulation  results  can  then 
be  used  to  develop  improved  closure  models  for  higher-order  correlations 
that  appear  in  the  models  for  the  Reynolds  stress  tensor. 

Within  this  framework,  examples  of  complex  aerodynamic  flow  fields 
including  a high-lift  system  composed  of  a multi-element  airfoil  configura- 
tion and  a scramjet/ramjet  engine  configuration  were  used  to  highlight  the 
need  for  well  chosen  unit  problems  that  would  isolate  the  key  dynamics 
associated  with  such  flows. 

Currently,  composite  methodologies  are  appearing  which  attempt  to 
utilize  both  the  LES  and  RANS-type  formulations.  Formal  and  complete 
methods  which  correctly  utilize  these  two  formulations  may  need  to  satisfy 
some  consistency  requirements  in  order  to  insure  that  both  the  LES  and 
RANS-type  formulation  will  yield  the  same  mean  flow  fields. 
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